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Abstract 

Based on the ab initio electronic structure calculations the picture of ferromagnetism in polimerized Ceo 
is proposed which seems to explain the whole set of controversial experimental data. We have demonstrated 
that, in contrast with cubic fullerene, in rhombohedral Ceo the segregation of iron atoms is energetically 
unfavorable which is a strong argument in favor of intrinsic character of carbon ferromagnetism which can 
be caused by vacancies with unpaired magnetic electrons. It is shown that: (i) energy formation of the 
vacancies in the rhombohedral phase of Ceo is essentially smaller than in the cubic phase, (ii) there is a 
strong ferromagnetic exchange interactions between carbon cages containing the vacancies, (iii) presence of 
iron impurities can diminish essentially the formation energy of intrinsic defects, and (iv) the fusion of the 
magnetic single vacancies into nonmagnetic bivacancies is energetically favorable. The latter can explain a 
fragility of the ferromagnetism. 

1 Introduction 

Carbon plays unique role in nature, from nucleosynthesis processes in stars to its crucial importance as a basic 
material substrate of life. It is not surprising therefore that all its properties is a subject of great interest. 
Magnetic properties of elemental carbon were studied in a context of geology and cosmology lj , biochemistry 
[2J, physics [3J, and material sciences [I]. The possibility of intrinsic long-range magnetic order in carbon is 
intriguing from both basic and applied scientific point of view since the properties of sp-electron magnetic 
materials can be essentially different from conventional itinerant-electron ferromagnets, in particular, it follows 
from general theoretical considerations that they may be perspective candidates to magnetic semiconductors 
with the Curie points Tq above room temperature [5]. 

The first observation [S] the ferromagnetism with Tq about 500K in polymerized fullerenes attracts much of 
attention and now this is the most studied magnetic carbon-based system. The experimental data are still very 
controversial and there is still no common opinion about intrinsic or extrinsic character of this ferromagnetism. 
The reported results were independently reproduced afterwards pj [5] and conditions of appearance of the 
ferromagnetism in polymerized fullerenes were studied in more detail EH EH EH EO] • It was found that 
the ferromagnetism arises only in a very narrow range of pressure and temperature in the vicinity of the point 
where fullerene cages are collapsed [2 EH] (for a general review of current status of the problem, see Ref. EI])- 
A typical criticism |14j which always occurs when discussing ferromagnetism in sp-electron systems with small 
magnetization (cf. the discussion for CaBg [5]) is a possible contamination of samples by iron and iron-based 
compounds due to use of metallic tools, however, a very careful experimental procedure avoiding this problem 
is described in Ref. E3] ■ The most serious criticism that have even resulted in the retraction of the first report 
ES] was presented in Ref. ES] where on the base of closeness of the Curie temperatures of ferromagnetic Ceo 
and cementite (FeaC) the effect was attributed to the presence of the latter phase. An idea of crucial role of 
Fe impurities and formation of ferromagnetic FesC phase can not explain the disappearance of ferromagnetism 
beyond a very narrow pressure and temperature interval of sample preparation [HI ESj since under high 
temperatures and pressures the cementite decays into kish graphite and pure iron, the latter having even higher 
Curie temperature than cementite. Also, it is very difficult to explain in these terms a strong dependence of the 
magnetization of photopolymerized fullerenes on the choice of oxygen-free or oxygen-rich formation conditions 
[TTj . The Curie temperature about 820K which was observed in one of the samples [8] does not coincide 
with those for cementite and a-iron. It is important to stress that, actually, only magnetic carbon systems 
synthesized from ethanol [17] can be considered as pure enough in a sense of iron impurities whereas natural 
materials such as coal and graphite always contain noticeable amount of iron and cannot be completely purified 
EH1 [201 [21]- Thus, nanotubes prepared from them are also contaminated by iron [52] and iron oxides [25] 
which are rather difficult to eliminate |24j . It is not surprizing therefore that sometimes an amount of iron in 
fullerene samples under discussion is as high as 400//gr/gr US] but sometimes it is negligible E3]. However, 
for the clearest samples in Ref. |13j ferromagnetism was not observed so, indeed, it is a temptation to connect 
the ferromagnetism of polymerized fullerenes simply with some iron-based magnetic phases, such as elemental 
iron, magnetite Fea04, or cementite FeaC E2]. However, this cannot explain the whole set of experimental data 
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available. First, magnetic characteristics of all these phases are well known and usually do not coincide with 
those of the "magnetic carbon" [TO], besides that, X-ray emission spectra of carbon in cementite are rather 
peculiar which allows us always to distinguish cementite from other carbon compounds [26) . The magnetic 
measurements was used already to separate intrinsic carbon magnetism from that of "parasytic" phases [1J (for 
polymerized fullerenes it was done in Ref. |10j). Also, even intentional introduction of iron to carbon systems 
is not sufficient to produce the ferromagnetism [571 HH] • 

In last few years new results appear about ferromagnetism in heavy-ion irradiated fullerene films [291 130] 
and proton irradiated graphite 31j. The x-ray magnetic circular dichroism measurements demonstrate mag- 
netic moments on carbon in proton irradiated graphite [35] . Another argument for intinsic character of carbon 
magnetism were discussed in Ref. |33j based on experimental data on C60H24. This provides some additional 
arguments in favor of possibility of intrinsic carbon ferromagnetism under certain conditions. To summarize 
the experimental data discussed above, one can conclude that there is a correlation between contamination by 
iron and ferromagnetism in the polymerized fullerenes [14] but it cannot be probably explained by the most 
straightforward way, as a result of formation of iron-based ferromagnetic phases. At least, some experimental 
data seem to be in disagreement with this view. Actually, the role of iron impurities in the problem of magnetic 
carbon may be more subtle and some more complicated opportunities can be considered, e.g., a formation of 
different molecular magnets based on charged fullerenes, with the Curie temperature Tc above room temper- 
ature (see Ref.[34 and references therein). In prinicple, this might explain a variety of Tc values in magnetic 
fullerenes as well as a fragility of their ferromagnetism. However, these systems [341 are periodic crystals and it 
is not known whether they can be ferromagnetic, with high enough Tc, under the disorder, or not. Also, they 
have specific Raman spectra different from those observed in the magnetic fullerenes. 

Another possible explanation is a catalytic activity of transition metals, in particular, iron, and their oxides, 
which are usually exploited for the nanotube gropwth |35[ I36[ I37j . One can assume that iron impurities can 
initiate a formation of various defects responcible for the ferromagnetism. It is the scenario which will be 
discussed further. 

Thus, the nature of magnetism in carbon based systems is still an open issue and thus theoretical inves- 
tigations can be very useful. Based on first-principle calculations of electron energy spectra and comparison 
with relevant experimental data it was shown [38j that the ferromagnetism cannot be an intrinsic property 
of ideal rhombohedral crystal lattice of the polymerized fullerenes. A role of doping by light elements in the 
appearance of ferromagnetism has been studied theoretically in Ref. [10] . A possible relation between vacancies 
and magnetism in carbon based systems was investigated intensively in last years. It is known |44( I45[ I46j 
that a single vacancy in graphene produces a local magnetic moment about 1.2 [Ib- At the same time, the 
vacancies in fullerenes were studied rather in a context of their role in structure stability of carbon cages [UJ [42] 
and information about their magnetic properties is rather poor. Magnetic interactions between vacancies in the 
polymerized fullerenes were studied in Ref. [43) but only for a one special kind of vacancies. The model proposed 
in Ref. cannot explain the ferromagnetism of proton- irradiated samples [2"9l [3Ti] , as well as conservation of 
ferromagnetism after the depolymerization [llj . At the same time, keeping in mind strong suspicions about an 
extrinsic character of "carbon" magnetism a theoretical clarification of issues on iron-carbon phases and metal 
impurities is necessary which was not done yet. 

Here we present the results of first principles calculations of magnetic and structural properties of various 
types of vacancies in fullerenes. It turns out that the recombination of magnetic single vacancies into non- 
magnetic bivacancies (which have the lower energy than a pair of the single vacancies) is the main limiting 
factor. This probably explains appearance of the ferromagnetism only in a narrow interval of pressure and 
temperature |13j which was not discusses in the earlier theoretical considerations. Based on our computational 
results we present also some arguments against decisive role of iron confirming that the ferromagnetism can be 
an intrinsic property of (defected) carbon systems. 



2 Computational method 

We used the SIESTA package for electronic structure calculations [J71 05] with the generalized gradient approx- 
imation for the density functional [3§J with energy cutoff 400 Ry, and fc-point 8x8x8 mesh in Monkhorst-Park 
scheme [50J . Coordinates of atoms and lattice parameters were taken for perfect (defect-free) cubic Cgo struc- 
ture from Ref. [5T], and for perfect tetragonal and rhombohedral phases from Ref. [S5]. Calculated electronic 
structure for both cubic and rhombohedral Cgo are in a good agreement with the experimental X-Ray spectra 
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3 Iron impurities on Ceo 



First, we investigated possible segregation of Fe impurities. To this aim, we have calculated solution (defect 
formation) energy for two Fe atoms situated in the same Cqq ball and its dependence on the distance between 
impurity for both cubic and rhombohedral phase, with the optimization of atomic positions in the buckyball. 
The lattice constants were kept fixed as it should be in the case of small impurity concentration. The defect 
formation energy was calculated as Ef orm = Ec B0 + n Fe — Ec 60 + nEp e , where Ec 60 + n Fe is the energy of the ball 
Ceo with adsorbed n atoms of iron, Ec 60 is the total energy of pure Ceo, and Ep e is the total energy of a-Fe per 
atom. Different geometric configurations for these calculations and the calculated energies are shown in Fig. [T] 
It is seen that for the cubic phase the single Fe atom per Ceo ball has the highest energy whereas for the pair 
of atoms the energy decreases with the distance decrease. This means that iron impurities in the cubic phase 
will segregate, indeed. However, for the rhombohedral phase the tendency is opposite. The formation energy in 
this case is the lowest for the single atom at the apex of the buckyball and for the pairs of Fe atoms it decreases 
with the distance increase. For the atomic configuration shown in Fig. [Tg the formation energy is smaller than 
for the single atom which is probably related with geometric distortions of the Cgo sphere, the iron atom being 
located in the void between neighboring buckyballs. We have calculated additionally the formation energy for 
the single Fe atom situated near the central part of rhombohedral carbon cage (Fig. [TJ) . This corresponds to the 
minimal energy. Thus, for the rhombohedral Cqo the segregation of Fe looks energetically unfavorable. Instead, 
at high concentration of impurities one can expect a formation of linear chains metal-buckyball-metal-... which 
were experimentally observed for the mixture of cobalt with Cgo [35] . 

4 Vacancies in Ceo 

The vacancy formation energy was calculated as Ef orm = E vac — Ec 60 (60 — n) / 60, where E vac is the total energy 
of the system with single vacancy per buckyball. As was shown in Ref. [I3J the vacancies can be responsible for 
ferromagnetism of fullerenes. The main problem is their too high energy formation, at least, in the cubic phase 
41 . However, the vacancy energy is sensitive to local geometric structure of fullerenes and, as we will show, 
there is an essential difference between cubic and rhombohedral phases. Figure [U shows a pictorial view of the 
cubic, tetragonal and rhombohedral phases. Whereas all atoms in the cubic phase are equivalent, already in the 
tetragonal one there are several different types of vacancy positions, due to formation of covalent bonds between 
the buckyballs and their distortions. We have performed calculations of the vacancy formation energy removing 
nonequivalent single atoms as shown in Figure [5] The calculated formation energy for the cubic phase is 5.58 
eV which is smaller than the value 7.2 eV for graphene [H]. In the tetragonal phase there are three types of the 
vacancy positions, the formation energy being decreased from the apex of the buckyball to its center. This is a 
consequence of the carbon-carbon bonds between the buckyballs in the central area which leads also to partial 
suppression of magnetic moments of these vacancies. The calculated formation energies for the tetragonal phase 
vary from 4.67 eV to 6.28 eV. Similar to Ref. [43] magnetic moments on vacancies are strongly dependent on 
geometric structure of fullcrcnc. 

The rhombohedral Cqq is characterized by a larger number of the C-C bonds and, as a result of polymer- 
ization, the ball is compressed by 8.4% along Z axis and expanded in the perpendicular directions and, thus, 
essentially different parts can be separated, the central one and two "capping" parts. The formation energy for 
vacancies in the capping parts (atoms 1 - 4 in Figure is about 5 eV which is close to corresponding values 
for other Ceo phases. However, for the central part (atoms 5-8 Figure ^jp) the vacancy formation energy is 
much smaller and for the atoms 5 and 7 it reaches 3.22 eV. That is why, probably, the polymerization at high 
pressure and high temperature is the way to magnetic carbon: it diminishes the vacancy formation energy and 
thus leads to a drastic increase in the vacancy concentration. 

The magnetism on the vacancies turns out to be rather fragile. The broken bonds required for magnetism 
are, at the same time, active centers for chemisorption. As a result, the magnetic moments can be destroyed 
by impurities. We have calculated the energy of oxygen chemisorption on a single vacancy. It is equal to -1.114 
eV for the most energetically favorable vacancy in the rhombohedral case. At the same time, for the vacancy in 
cubic phase this energy is positive and equals 0.716 eV. These data can explain suppression of ferromagnetism 
at the photopolimerization of Cgo under oxygen-rich conditions |llj . 

Figure [3] presents electron densities of states for all studied phases of Cgo with and without vacancies. At a 
large energy scale, vacancies modify the electronic structure rather weakly so it would be difficult to see their 
effects in, say, X-ray spectra [31] . At the same time, they change essentially the electronic structure in the close 
vicinity of the Fermi level producing conduction electrons in the former energy gap. Note that experimentally 
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Figure 1: Optimized geometrical structure of iron atoms adsorbed on Ceo in cubic and rhombohedral phases. Numbers 
near the pictures are defect formation energies per iron atom, in eV. 
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Figure 2: Top (left) and side (right) view of cubic (a), tetragonal (b) and rhombohedral (c) fullerenes. Red (dark) 
spheres show inequivalent atoms which can be removed. The values of magnetic moments per whole buckyball are shown 
in parentheses (in /zb). 
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Figure 3: Total densities of states for (a) cubic, (b) tetragonal, (c) rombohedral fullerens. Solid red and dashed green 
lines correspond to the system without and with vacancy, respectively. 
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Figure 4: Distribution of magnetic moments in jib on Ceo sphere with single vacancy in rhombohedral phase. 

magnetic samples turned out to be conducting |12j . These charge carriers may be responsible for ferromagnetism 
appearance in a narrow defect band [5] . 

To estimate the exchange interaction energy responsible for possible magnetic ordering we have performed a 
total energy calculation for a supercell containing two Ceo balls with single vacancy per ball, with parallel (FM) 
and antiparallel (AFM) spin orientations. The calculations have been done for both cubic and rhombohedral 
phase, in the latest case the vacancies with the lowest formation energy being considered. For both systems the 
atomic relaxation has been taken into account. We have found that the difference between total energy in FM 
and AFM configurations is 0.41 eV for the rhombohedral phase and 0.53 eV for the cubic one. These values 
are close to those calculated for the hydrogen-doped Ceo [10] ■ This energy difference is very high which can 
be explained by a peculiar distribution of the electron spin density. Whereas in open structures like graphene 
the magnetic moments from vacancies and impurities are rather localized, with a fast decay of the spin density 
[22 US] , in closed structures such as Ceo the magnetic moments distribute over the whole buckyball (see Figure 
IU cf. Ref. [40] for the hydrogen-doped fullerenes). Of course, the value of exchange interactions between distant 
balls with the vacancies will be weaker than that between the neighboring ones so our calculation gives just an 
upper limit for the value of magnetic interactions. Various concentration of fullerenes with vacancies will result 
in various distances between them and, thus, with various exchange parameters and Curie temperatures which 
seems to be in agreement with experimental data showing a large variance of Tc in polymerized fullerenes, from 
550 K to 800 K. 

Let us consider now an effect of iron impurities on the vacancy formation in fullerenes. By analogy with the 
process of ripping of few layer graphene by iron nanoparticles [55] the formation energy will be calculated by 
the expression Ef orrn = Ec 60 Fe — Ec 59 Fe — Ec 60 59/60. As the simplest model we consider the system C^Fe 
where carbon atom closest to the iron impurity is removed. We start with the cubic phase of fullerene. The 
computational results are presented in FigHK One can see that in the presence of Fe the vacancy formation 
energy is drastically decreased. Since the iron atoms in cubic phase tend to segregate as was discussed above 
this caluclation deals with rather hypopthetic situation than real one. However, a similar tendency holds also 
for the polymerized fullerenes. As shown in Fig[5jD-d for all positions of Fe under consideration the vacancy 
formation energy is essentially lower than without iron. Because single iron atoms without magnetic carbon 
impurities in neighbors stay paramagnetic centers |28j . Thus, contamination by iron can initiate a formation 
of defects responsible for ferromagnetism. The effect of Fe on the vacancy formation cosidered here may be 
also interesting in light of simulation of metal-fullcrcnc systems for potential application in hydrogen storage 



Based on the computational results one can conclude that vacancies can be a cause for the ferromagnetism 



57, 56J. 
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Figure 5: Optimized structures of fullerenes with iron impurity and carbon vacancy nearby for cubic phase (a) and 
nonequivalent positions in rhombohedral (b)-(d). The numbers correspond to formation energies of the vacancy in 
electronvolts with (without) Fe impurity. 
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Figure 6: Optimized geometrical structures of mono-, bi-, tri- and quadrovacancies in rhombohedral Ceo. 




12 3 4 

Number of removed atoms 

Figure 7: Energy formation for multivacancies per removed atom as a function of the number of removed atoms in 
central part of rhombohedral Ceo (solid red line), on top of the rhombohedral Ceo (dashed green line), and in cubic Ceo 
(dotted blue line). 
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of polymerized fullerenes. To understand, in a framework of this picture, why the ferromagnetism is so fragile 
and why it disappears at the temperature of pressure increase the energetics of multivacancy formation has been 
considered. We have calculated the energies and magnetic moments for bi-, tri- and quadrovacancies for both 
cubic and rhombohedral Ceo (see Fig. [SJ. The computational results presented in Fig. [7] demonstrate that the 
formation energies per removed atom for the bivacancies are smaller than for the single vacancy, similar to the 
previous results for graphene [SI]. Therefore, the process of vacancy fusion into the bivacancy is energetically 
profitable. The bivacancy is always nonmagnetic since there is no unpaired electrons in the situation of two 
broken bonds. This explains a very narrow temperature and pressure interval where the ferromagnetism in 
polymerised Cgo arises: further increase of the pressure or temperature stimulates the process of the bivacancy 
formation. Figure [7] shows that the appearance of the multivacancies is the most probable in the central part 
of rhombohedral Cgo • This is exactly the place where the cage collapse happens as is confirmed by observation 
of curved fragments of fullerenes in the "hard carbon" after the cage collapse [HJ [59] . 

5 Conclusion 

Thus, the explanation of appearance and disappearance of ferromagnetism in polymerized fullerenes in terms 
of vacancies seems to be consistent with the whole set of experimental data. According to this picture, the 
ferromagnetism arises under rather restrictive conditions where the monovacancies are formed but their fusion 
into nonmagnetic bivacancies can be prevented. Oxygen chemisorption will also suppress the ferromagnetism due 
to formation of chemical bonds involving unpaired electrons associated with the vacancies. We have discussed 
also mechanisms of the cage collapse via fromation of quadrocvacancy in the central part of buckyballs. We have 
demonstrated that a segregation of iron impurities in the rhombohedral fullerene is energetically unprofitable, 
in contrast with the cubic one. Importantly, a contamination by iron can decrease the vacancy formation energy 
in fullerenes. 

The work is financially supported by Stichting voor Fundamenteel Onderzoek der Materie (FOM), the 
Netherlands. 

References 

[I] J. D. M. Coey, M. Venkatesan, C. B. Fitzgerald, A. p. Douvalis, I. S. Sanders, Nature 420 (2002) 156-159. 
[2] E. Coronado, J. R. Galan-Mascaros, C. J. Gomes-Garsfa, V. Laukhin, Nature 408, (2000) 447-449. 

[3] T. W. Odom, J. L. Huang, C. L. Cheung, C. M.Lieber, Science 290, (2000) 1549-1552. 

[4] L. Krusin-Elbaum, D. M. Newns, H. Zeng, V. Derycke, J. Z. Sun, R. Sandstrom, Nature 431, (2004) 672-676. 

[5] D. M. Edwards, M. I. Katsnelson, J. Phys.: Condons. Matter 18, (2006) 7209-26. 

[6] T. L. Makarova, B.Sundqvist, P. Esquinazi, Y. Kopelevich, P. Scharff, V. A. Davydov, L. S. Kashevarova, 
A. V. Rakhmanina, Nature 413 (2001) 716-9. 

[7] R. A. Wood, M. H. Lewis, M. R. Lees, S. M. Bennington, M. G. Cain, N. Kitamura, J Phys: Condens 
Matter 14, (2002), L385-392. 

[8] V. N. Narozhnyi, K.-H. Miiller, D. Eckert, A. Teresak, L. Duunsch, V. A. Davydov, L. S. Kashevarova, A. 
V. Rakhmanina, Physica B 329-233, (2003) 1217-1218. 

[9] T. L. Makarova, B. Sundqvist, Y. Kopelevich, Synth Metals 137, (2003), 1335-1337. 

[10] T. L. Makarova, I. B. Zakharova, Fullerenes, Nanotubes and Carbon Nanostructures 16, (2008) 567-573. 

[II] T. L. Makarova, K.-H. Han, P. Esquinazi, R. R. da Silva, Y. Kopelevich, I. B. Zakharova, B. Sundqvist, 
Carbon 41, (2003) 1575-1584. 

[12] T. L. Makarova, B. Sundqvist, P. Schaff, M. E. Gaevski, E. Olsson, V. A. Davydov, A. V. Rakhmanina, L. 
S. Kashevarova, Carbon 39, (2001) 2203-2209. 

[13] K.-H. Han, A. Talyzin, A. Dzwilewski, T. L. Makarova, R. Hone, P. Esquinazi, D. Spemann, L. S. Dubrovin- 
sky, Phys. Rev. B 72 (2005) 224424. 



10 



[14] A. V. Talyzin, A. Dzwilewski, J. Nanoscience Nanotechnol 7, (2007) 1151-1161. 

[15] A. V. Talyzin, A. Dzwilewski, L. Dubrovinsky, A. Setzcr, P. Esquinazi, Eur. Phys. J. B 55, (2007) 57-62. 

[16] T. L. Makarova, B.Sundqvist, P. Esquinazi, Y. Kopelevich, P. Scharff, V. A. Davydov, L. S. Kashevarova, 
A. V. Rakhmanina, Nature 440 (2006) 707. 

[17] N. Parkansky, B. Altcrkop, R. L. Boxman, G. Leitus, O. Berkh, Z. Barkay Z, Yu,. Rosenberg, N. Eliaz, 
Carbon 46 (2008) 215-219. 

[18] K. K. Bcrtinc, E. D. Goldberg, Science 173, (1971) 233-235. 

[19] W. Primak, L. H. Fuchs, Phys. Rev. 95, (1954) 22-31. 

[20] Y. Kaburagi, Y. Hishiyama, Carbon 36, (1998) 1671-1676. 

[21] P. Esquinazi, A. Setzer, R. Holme R, C. Semmclhack, Y. Kopelevich, D. Spemann, T. Butz, B. Kohlistrunk, 
M. Losche, Phys. Rev. B 66, (2002) 024429. 

[22] K. A. Williams, M. Tachibana, J. L. Allen, L. Grigorian, S.-C. Cheng, S. L. Fang, G. U. Sumanasekera, A. 
L. Lopcr, J. H. Williams, P. C. Eklund, Chcm. Phys. Lett. 310, (1999) 31-37. 

[23] J. Kruusma, N. Mould, K. Jurkschat, A. Crossley, C. E. Banks, Elcctrochem. Commun. 9,(2007) 2330-2333. 

[24] K. Jurkschat, X. Ji, A. Crossley, R. G. Compton, C. E. Banks, Analyst 132, (2002) 21. 

[25] R. Hohnc, P. Esquinazi, Adv. Mater. 14, (2002) 753-756. 

[26] E. Z. Kurmacv, S. N. Shamin, K. M. Kolobova, S. V. Shulepov, Carbon 24, (1985) 249-253. 

[27] J. Barzola-Quiquia, R. Hohne, M. Rothermel, A. Setzer, P. Esquinazi, V. Hccra, Eur. Phys. J. B 1, (2008) 
127-130. 

[28] R. Siclcmann, Y. Kobayashi, Y. Yoshida, H. P. Gunnlaugsson, G. Weyer, Phys. Rev. Lett 101, (2008) 
137206. 

[29] A. Kumar, D. K. Avasthi, J. C. Pivin, A. Tripathi, F. Singh, Phys. Rev. B 74, (2006) 153409. 

[30] A. Kumar, D. K. Avasthi, J. C. Pivin, R. M. Papaleo, A. Tripathi, F.Singh, I. Sulania, J. Nanosci Nan- 
otechnol 7, (2007) 2201-2205. 

[31] P. Esquinazi, D. Spemann, R. Hohnc, A. Setzer, K.-H. Han, T. Butz, Phys. Rev. Lett. 91, (2003) 227201. 

[32] H. Ohldag, T. Tyliszczak, R. Hohnc, D. Spemann, P. Esquinazi, M. Ungureanu, T. Butz, Phys. Rev. Lett. 
98, (2007) 187204. 

[33] V. E. Antonov, I. O. Bashkin, S. S. Khasanov, A. P. Moravsky, Yu. G. Morozov, Yu. M. Shulga, Yu. A. 
Osipyan, E. G. Ponyatovsky, J. Alloys and Compounds 330-333 (2002) 365-368. 

[34] D. V. Konarcv, S. S. Khasanov, A. Otsuka, G. Saito, R. N. Lyubovskaya, J. Am. Chcm. Soc. 128, (2006) 
9292-9293. 

[35] D. Takagi, Y. Homma, H. Hibino, S. Suzuki, Y. Kobayashi, Nano Lett. 6, (2006) 2639-2641. 

[36] Zhou W, Han Z, Wang J, Zhang Y, Jin Z, Sun X, Y. Zhang, C. Yan, Y. Li, Nano Lett 6, (2006) 2987-2990. 

[37] O. V. Yazyev, A. Pasquarello, Phys. Rev. Lett. 100, (2008) 156102. 

[38] D. W. Boukhvalov DW et al., Phys. Rev. B 69, (2004) 115425. 

[39] V. Lavrentiev, H. Abe, S. Naramoto, S. Sakai, K. Narumi, Chcm. Phys. Lett. 424, (2006):101-104. 

[40] O. E. Kvyatkovskii, I. B. Zakharova, A. L. Shclankov, T. L. Makarova, Phys. Rev. B 72, (2005) 214426. 

[41] Y.-H. Kim, J. Choi, K. J. Chang, Phys. Rev. B 68 (2003) 125420. 

[42] Y.-H. Hu, E. Ruckcnstcin, J. Chcm. Phys. 119, (2003) 10073. 



11 



[43] A. N. Androtis, M. Menon, R. M. Sheez, Phys. Rev. Lett. 90, (2003) 026801. 

[44] P. O. Lehtinen, A. S. Foster, Y. Ma, A. V. Krasheninnikov, R. M. Nicmicn, Phys. Rev. Lett. 93, (2004) 
187202. 

[45] E. J. Duplock, M. Scheffler, P. J. D. Lindan, Phys. Rev. Lett. 92, (2004) 225502. 
[46] O. V. Yazyev, L. Helm, Phys. Rev. B 75, (2007) 125408. 

[47] E. Artacho, J. D. Gale, A. Garsia, J. Junquera, R. M. Martin, P. Orejon, D. Sanchez-Portal, J. M. Soler, 
SIESTA, Version 1.3 2004. 

[48] J. M. Soler, E. Artacho, J. D. Gale, A. Garsia, J. Junquera, P. Orejon, D. Sanchez-Portal, J. Phys.: 
Condons. Matter. 14, (2002) 2745-2781. 

[49] J. P. Perdcw, K. Burke, M. Ernzcrhof, Phys. Rev. Lett. 77, (1976) 3865-3868. 

[50] H. J. Monkhorst, J. D. Park, Phys. Rev. B 13, (1976) 5188-5192. 

[51] S. Satapathy, V. P. Antropov, O. K. Andersen, O. Jepsen, O. Gunnarsson, A. I. Liechtenstein, Phys. Rev. 
B 46, (1992) 1773-1793. 

[52] M. Nuhez-Regueiro, L. Marques, J.-L.Hondeau, O. Bethoux, M. Perroux, Phys. Rev. Lett 74, (1995) 
278-281. 

[53] L. Pisani, B. Montanari, N. M. Harrison, New J. Phys. 10, (2008) 033002. 
[54] J. M. Carlson, M. Scheffler, Phys Rev Lett 96,(2006) 046806. 

[55] S. S. Datta, D. R. Strachan, S. M. Khamis, A. T. C. Johnson, Nano Lett. 8, (2008) 1912-1915. 
[56] S. Yang, M. Yoon, E. Wang, Z. Zhang, J. Chem. Phys. 129, (2008) 134707. 

[57] Y. Zhao, Y.-H. Kim, A. C. Dillon, M. J. Hcbcn, S. B. Zhang, Phys. Rev. Lett. 94, (2005) 155504. 

[58] E. V. Tat'yanin, A. G. Lyapin, V. V. Mukhamadiarov, V. V. Brazhkin, A. L. Valilicv, J. Phys.: Condens. 
Matter 17, (2005) 249-257. 

[59] V. V. Brazhkin, V. L. Solozhcnko, V. I. Bugakov, S. N. Dub, O. O.Kurakevych, M. V. Kondrin, A. G. 
Lyapin, J. Phys.: Condens. Matter 19, (2007) 236209. 



12 



